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Abstract 

Motivated by the theoretical possibility of floating orbits and the potential to contribute extra 
constraints on alternative theories, in this paper we derive the self-force equation for a small 
compact object moving on an accelerated world line in a background spacetime which is a solution 
of the coupled gravitational and scalar field equations of scalar-tensor theory. In the Einstein frame, 
the coupled field equations governing the perturbations sourced by the particle share the same form 
as the field equations for perturbations of a scalarvac spacetime in general relativity, with both 
falling under the general class of hyperbolic field equations studied in |l|]. Here, we solve the field 
equations formally in terms of retarded Green functions, which have explicit representations as 
Hadamard forms in the neighbourhood of the world line. Using a quasi-local expansion of the 
Hadamard form, we derive the regular solutions in Fermi normal coordinates according to the 
Detweiler-Whiting prescription. To compute the equation of motion, we parameterize the world 
line by the particle’s mass and “charge”, which we define in terms of the original Jordan frame 
mass, its derivative, and the parameter which translates the proper time in the Jordan frame to the 
Einstein frame. These parameters depend on the value of the background scalar field and its self- 
field corrections. The equation of motion which follows from the regular fields strongly resembles 
the equation for the self-force acting on a charged, massive particle in a scalarvac geometry of 
general relativity. Unlike the scalar vacuum scenario, the “charge” parameter in the scalar-tensor 
self-force equation is time variable and leading to additional local and tail terms. We also provide 
evolution equations for the world line parameters under the influence of the self-fields. 

I. INTRODUCTION 

A. Motivation for alternative theories 

The general theory of relativity (GR) is widely considered to be the most successful 
classical theory of gravitational phenomena. Its predictions have passed many stringent 
tests such as the light deflection around the sun during a solar eclipse, the Shapiro time 
delay, the perihelion advance of Mercury, and the Nordtvedt effect in lunar motion j^. 
Despite its great success in predicting deviations from Newtonian gravity in the solar system 
where deviations from flat spacetime are weak and non-relativistic {v ^ 10“^c, and U/A ~ 
10“®), it has received little direct experimental verihcation in the strong-held regime, which 
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leaves room for alternative theories of gravity that reduce to general relativity in the weak- 
held limit. Additional evidence supporting general relativity has come from observations 
of binary pulsar systems which measure several Keplerian and post-Keplerian parameters 
including the decay of the orbital period of the binary due to binding energy lost in the 
form of gravitational waves. The parameters of some alternative theories which predict 
dipole radiation in addition to the quadrupole radiation predicted by general relativity have 
been constrained by observations of the neutron-star-white-dwarf systems J11416545 and 
J1738-I-0333 to an even higher degree than the bounds from the Cassini test in the solar 


system 




Research in alternative theories of gravity has been traditionally fueled in large part by 
the long standing theoretical effort to formulate a quantum theory of gravity, which promises 
to shed light on big mysteries such as the origin of the universe and the singularities hidden 
deep within the event horizons of black holes. A unihed theory of gravity and quantum 
mechanics is necessary because the general relativistic picture of spacetime as a Riemannian 
manifold breaks down at the length scale £p ~ 10“^^ cm where quantum fluctuations of 
the gravitational held become important. However, deviations from classical relativity may 
be exhibited in the strong-held regime even at macroscopic length scales. For instance, 
deviations from Einsteinian gravity at large scales may account for the dark energy content in 
the universe and active research includes using observational data to constrain models which 

n 

parameterize such deviations in structure formation models [5|. Furthermore, models with 
additional scalar helds play a particularly important role in modern theories of cosmology 
by providing the muscle for the inhationary epoch. 

Among the many alternative theories of gravity, a particularly well studied and well 
established theory is scalar-tensor (ST) gravity, wherein the agents of the gravitational force 
are a scalar and a tensor held. ST theories originate from the work of Jordan ml, Fierz 

n f] 

[7[, Brans and Dicke [8[, who sought a way to incorporate Mach’s principle into a covariant 
theory of gravity. They found that a dynamical gravitational constant built from a scalar 
held 0 ~ G~^ acted as a way to inhuence local physics through large scale phenomena, 
while upholding the weak equivalence principle and maintaining general covariance. The 
held equations of the theory they developed proved to be equivalent to a sector of the 
dimensionally reduced hve-dimensional held theory designed by Kaluza {q] as an attempt 
to unify gravity and electromagnetism. The emergence of a scalar held(s) is now known 
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to be a generic property of dimensionally reduced higher-dimensional models such as string 
theory, making the study of gravitational theories with additional scalar helds a very active 
research helc. 


Hawking 


34l | has shown that stationary solutions of Brans-Dicke type scalar tensor held 


theories must also satisfy the Einstein held equations and that the general stationary solu¬ 


tions are the Kerr family of so 
on this, Sotiriou and Faraoni 


utions coupled to a constant background scalar held. Building 


35| extended Hawking’s proof to include more general ST the¬ 


ories including f{R) theories and relaxed the need for a symmetric collapse. However, the 
proof fails to hold for certain non-stationary cosmological scalar helds, such as those which 
asymptote to solutions with linear time dependence at large radii 3^ , leading to observable 


deviations from general relativity in the orbital decay of super-massive black hole binaries 


37l |. Although constraints on such theories do exist, for generality we choose to leave the 


spacetime dependence of the background scalar held unspecihed and we make no restric¬ 
tions on the behaviour of the scalar potential acting as the cosmological function other than 
convexity. We then consider a stationary background scalar held conhguration as a special 


case. 


B. Compact binary problem 

With the promise of gravitational wave observatories directly probing general relativity 
in the strong-held regime through gravitational wave (GW) signals, researchers have begun 
to investigate the dynamics of compact binary systems like binary black holes and neutron 
stars within various alternative theories of gravity. 

Scalar-tensor theories beneht from a well-posed Cauchy formulation, enabling numerical 
relativity simulations to explore strong-held gravitational ehects [lO, 11|. Although post- 
Newtonian and linearized theory results indicate that binary black holes in ST theory are 


indistinguishable from their GR counterparts 


0,0, £i 


just like t 


in the general proofs of Hawking and Faraoni and Sotiriou 3^, 1 m 1, if the scalar held is 


re isolated black holes 


made time dependent or given some inhomogeneity through an external mechanism like a 
potential, the binary also emits dipolar radiation while accreting the scalar held. Healy and 
collaborators [0] used a Mexican hat type potential to induce signihcant non-uniformity in 


the scalar held to investigate the gravitational and scalar radiation prohles from merging 
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binary black holes in ST theory. Binary neutron star systems, on the other hand, are 
more sensitive to the presence of the scalar held and display new phenomena that are 
qualitatively different from predictions of GR. Isolated neutron stars undergo a phenomenon 


known as spontaneous scalarization 


12j, where they develop a sudden scalar charge when 


they previously had none. Spontaneous scalarization occurs when a non-zero value of the 
scalar held inside the star becomes energetically favorable over the zero held conhguration 
for certain values of the ST parameter u. In the non-linear binary merger regime, 111 
have shown that spontaneous scalarization occurs in much the same way when the merging 
neutron stars are sufficiently close. The same group has shown that the gravitational wave 
signals from binary neutron stars may provide constraints on the parameters of ST gravity 
if one of the companion stars becomes scalarized prior to the time when the gravitational 
waves enter the detectors sensitivity band H- 

When the binary objects are of similar size and in the early stages of the inspiral, the 
helds are relatively weak and the compact bodies move with small orbital velocity relative 
to the speed of light. Typical examples of these systems include binary neutron stars and 
stellar mass black holes which are several hundred Schwarschild radii apart. In this regime, 
post-Newtonian theory provides the framework for computing the orbital motion and grav¬ 


itational wave signals. Starting with the work of Wagoner 


16| . researchers have used the 


ling bil iary 
3, 17-20|. 


post-Newtonian approximation to derive equations of motion for slowly inspiral 
systems emitting gravitational waves using the held equations of ST theory 
These works express the period decay formula in terms of Brans-Dicke parameters and use 
it to place constraints on the theory through observations in the strong held regime. The 
post-Newtonian waveforms in ST theory have been shown to reduce to the general rela¬ 
tivity waveforms for binary black holes to second post-Newtonian order in the radiation 
sector 24, 251. Recent interest in cosmological theories with light scalar helds 0, Q and 
the existence of boating orbits 0]; motivated Alsing and collaborators 231 to explore the 
post-Newtonian regime of compact binaries with a massive scalar held leading to an explicit 
expression for the decay of the binary’s orbital period in massive Brans-Dicke theory. 

For compact binaries with largely disparate mass parameters, such as a stellar mass 
(m ~ Mq) object orbiting a super-massive black hole in the centre of a galaxy Msmbh ~ 
IO^Mq, the binary’s orbital velocity may enter the highly relativistic regime where post- 
Newtonian methods no longer work. Instead, one must solve the perturbed Einstein held 
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equations without any recourse to slow-motion expansions. Yunes and collaborators 261 
used a test particle approximation to calculate the gravitational wave energy flux coming 
from an EMRI system consisting of a small body orbiting a Kerr black hole to examine 
whether low-frequency gravitational wave detections may provide constraints on the Brans- 
Dicke parameters. They showed that for massless scalar fields, the constraints are notably 
worse than those provided by solar-system tests (dephasing of GWs is weakly effected by 
changes in the parameters), but that massive scalar fields can provide significant constraints 


27|. 


The main motivation for this work is to provide an additional framework to constrain 
these theories in the EMRI regime via the effects of the self-force. The prospect is quite 
good for nontrivial configurations of the scalar field, due to the presence of new local terms 
and of several new non-diagonal tail terms that appear in the equations of motion. We will 
also show that constraints may be promising even when the scalar field configuration is the 
trivial one, as in the black hole scenario, due to the existence of a scalar component to the 
self-force in addition to the gravitational component. The framework will also serve as a 
means to dynamically study the floating orbits which arise in motion of particles around 
spinning black holes in scalar-tensor theories. 

Floating orbits are a byproduct of the superradiant instability displayed by rotating black 
holes with massive external scalar fields where the net energy flux at infinity is entirely due to 
the rotational energy extracted from the spinning black hole, halting the inspiral of the small 
body; i.e., E implying Eorh = 0. In this scenario, scalar perturbations in the 
superradiant regime with frequencies close to the mass of the scalar held create resonances 
which lead to large amounts of negative scalar flux down the horizon balancing the flux of 
gravitational radiation at null infinity. The mass of the scalar held is necessary to act as a 
reflecting wall which traps the radiation in order to create quasi-bound states of the held. 
The arguments for the existence of floating orbits rely heavily on an adiabatic approximation 
in which the orbit of the small body is modeled as a sequence of quasi-circular geodesic orbits 
which get progressively smaller in radius as orbital binding energy is released in the form 
of gravitational waves. Alternatively, the power radiated by the particle causing its inspiral 
may be written in terms of the work done by the dissipative piece of the radiation reaction 
self-force per unit time. The self-force also contains a conservative component acting on 
orbital timescales which may effect the stability of the floating orbits. A full calculation of 
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the self-force in the EMRI problem within ST theories is necessary to study how floating 
orbits arise dynamically and to assess their stability, which is a signihcant motivator for the 
present work in addition to the general study of motion in alternative theories. 

The leading order equations of motion for small compact objects in ST theories were given 
by Gralla 


, who used the Bianchi identity and various universal scaling relations for the 


helds in the region outside a small body of mass M and charge Q to map the near-zone 
Coulomb-type behaviour of the perturbed helds to a set of effective distributional sources. 
With this he was able to rigorously show that the linearized held equations for a small 
extended body in Einstein frame ST theory reduced to the point particle form 


SGa /3 — S7l6Tai3 = Sn J M{T)UaUp5i{x,z{T))dT, 

= —StT [ Q{t) S 4 {x, z{T))dT, 


(la) 

(lb) 


where 6 is the perturbation operator, Ga/s is the Einstein tensor. Tap is the stress-energy 
tensor of the bulk scalar held, and □ := is the covariant wave operator in the back¬ 

ground spacetime. Gralla also noted that the charge Q of the small body is not constrained 
by any evolution equation inherent to ST theory; it must be prescribed from the model de¬ 
scribing the internal structure of the object. He went on to show that the mass and charge 
of the small body obey surface integral relations similar to the ADM expressions 29|]. The 
foundational treatment given by Gralla set the stage for the simultaneous works of Zimmer- 


Q 


man and Poisson 111 and Linz and collaborators 


32| . who overcame the technical challenge 


of solving the coupled held equations and derived equations of motion for the coupled self¬ 
force and regularization parameters in the scalarvac and electrovac spacetimes of GR. These 
derivations rely on the Detweiler-Whiting axioms to identify and remove the divergence in¬ 
troduced by the point-particle. In vacuum GR, self-consistent self-force derivations using 


30|, 


hnite sized bodies and matched-asymptotic expansions performed by Gralla and Wald 
and also Pound [^, have reproduced the results of axiomatic point-particle derivations. We 


expect that our axiomatic derivation will also agree with future non-vacuum self-consistent 
computations. 

Here we use the technology developed in the aforementioned papers to compute the self¬ 
force in generic massive scalar-tensor theories and we hnd that the resulting equations of 
motion are substantially diherent from those in vacuum GR due to the presence of the 
additional scalar held. 
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The main results of this work are the equations of motion for a small body due to the 
coupled self-force in a general ST theory in the Einstein frame. We hrst give the equation of 
motion for a background scalar field with arbitrary time dependence and a general potential 
function. For this we hnd that the particle’s charge can be expressed in terms of the ST 
parameters in such a way as to resemble the scalarvac self-force plus additional corrections. 
We then analyze the asymptotically flat black hole scenario where the scalar field is constant 
and the potential and its hrst derivative are zero. In all cases we hnd that the equation of 
motion can be decomposed into the form 


Ma'‘ = FS + F" + F",„ 


where Fq is the force resulting from the gradient of the background potential, is the 
local contribution to the self-force which is built from background quantities evaluated on 
the world line, and is the non-local contribution to the self-force which takes the form 
of a time integral over the particle’s past history. The absence of the coupling in the 
asymptotically hat, constant scalar held black hole case leads to a self-force which is simply 
the sum of the gravitational and scalar self-forces. 

The paper begins in Sec Hi] with a derivation of the perturbed held equations in the 


Einstein frame for the general ST theo^ of a point particle. In Sec IIII A1 we review the 
condensed index notation presented in [ij and introduce several quantities to facilitate the 


local expansion of the perturbed helds. Next, in Sec IIII B1 we derive an expression for the 
regular held and compute its gradient as required for the self-force. Finally, we derive the 
self-force equation of motion for a point-particle in ST theory (Sec. IIVP and specialize it to 
a black hole in an asymptotically hat spacetime (Sec. jVj). 


II. THE SCALAR-TENSOR THEORY OF A SINGLE SCALAR FIELD 


The quantity of information on ST theories is substantial, but here we limit ourselves to 


introducing elementary quantities and concepts of ST theory needed for t 
self-force derivation. A detailed review of the theory may be found in 


le purposes of the 


The development here is to start with the action in the Jordan frame and then perform 
a conformal transformation to the Einstein frame where the held operator is the Einstein 
tensor and the scalar held contributes only as matter content through the bulk stress-energy 











tensor. We then vary the action with respect to the Einstein frame helds to obtain the held 
eqnations. This allows ns to express the pertnrbed held eqnations in manifestly hyperbolic 
form after a snitable gange transformation. From there we can express the pertnrbed held 
eqnations in a form where we can directly apply the mnlti-held methods of [l|. 

The action for a generic scalar-tensor theory in the Jordan frame has the form 


^.7 = 


Ibvr 


a{^)R - - 2c(0)) + 5 m(Tm, ^1, 


( 2 ) 


where Tm collects the helds responsible for the matter content of the theory, a, b and c are 
held-dependent ST parameters, and the overbar indicates that the qnantity is taken with 
respect to the Jordan frame. Withont loss of generality, we can redehne the scalar held 
a(0) 0. After the held redehnition we see that the most general ST theory contains two 

free fnnctions: a conpling fnnction (typically written as u}{4>)) which precedes the kinetic 
term and a cosmological fnnction which enters withont any explicit powers of the scalar 
held. The conpling fnnction is responsible for the spontaneons scalarization phenomenon 
mentioned in the introdnction where the binary develops a sndden scalar charge 2^. The 
ehects of the cosmological fnnction inclnde providing the scalar held with mass and playing 
the role of the cosmological constant. The potential also contribntes corrections to the mass 
evolntion of scalar particles in cnrved spacetime as we shall see. 

In the EMRI scenario of vacnnm general relativity, the matter content is a small body 
which is approximately described by a strnctnreless point-particle of mass m. In the ST 
theory of material objects, however, the inertial mass and strnctnral properties of the small 
body will generally be inhnenced by the scalar held dne to the variability of Newton’s grav¬ 
itational constant. This creates a need to introdnce a new parameter which parameterizes 
the sensitivity of the body’s binding energy on the backgronnd scalar held. The sensitivity 
of nentron stars will depend on their radins and eqnation of state bnt black holes all have the 
same sensitivity. The dependence on the scalar held is incorporated into the point-particle 
model by allowing the mass to vary with the scalar held 
the theory is thns given by 


. The point-particle action for 


^M = 


m{(p) dr 


(3) 


in the Jordan frame. Althongh the Jordan frame benehts from explicit conservation of 
stress-energy VaTj^ = 0, the held eqnations in the Jordan frame do not lend themselves 


to a direct application of the methods of Zimmerman and Poisson |l|. For this, we mnst 
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transform the fields into the Einstein frame. Performing the conformal transformation [1 


9fj,y 

X(0) 


f (3 /aX0)\ 

J UU(0); 


Mx) = a 


Fix) 


c(0) 

a2(0)’ 


we are led to the Einstein frame action 


2a{(p) j 


d(j), 


(4a) 

(4b) 

(4c) 

(4d) 


/ (i? - 9^^^,X^.X - 2F(x)) dV- j A{x)m{x)dT, (5) 

where we have introdnced the covariant volnme element in the Einstein frame dV : = 
y/—gd‘^x. A variation with respect to the metric yields the held eqnations in the Einstein 
frame 

G.„ = 8>r (r^y + T") , (6) 

where 

S (VaX'Wx + 2F)j (7) 

is the stress-energy of the bnlk scalar held and 

= J A(x)m{x)Uf,u^S4ix,z)dT (8) 

is the stress-energy of the point particle. Using the Einstein held eqnation, we hnd that the 
Ricci tensor in terms the matter held reads 


Rais = ^aX^fSX + FgajS- (9) 

Varying the action with respect to the scalar held prodnces the scalar wave eqnation 

□x - F'{x) = j ^^^^^ <^ 4 ( 3 ^, z)dT (10) 

governing the evolntion of y, where □ := g^'^'W is the covariant wave operator and 
F'ix) •= dF/dx- The presence of the point-particle creates pertnrbations of the helds 
aronnd their backgronnd valnes. If we denote the backgronnd helds (fnll helds taken at 
m = 0) by x(0) := $ and g^y{R) := g^p, the pertnrbed helds are given by 


/:=X-$, 

hfip ■ g^u 9^J,l2■ 
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(lla) 

(llb) 







In what follows we always work with either the metric perturbation or the background 
metric, so we drop the tilde notation and refer to the background metric simply as gafi- It 
is convenient to work with the trace-reversed metric perturbation 

where h = The beneht of the trace-reversed metric is that it has zero divergence 

in the Lorenz gauge. Here we require the metric perturbation to satisfy the one-parameter 
family of gauge conditions 

= 2XfV^^ 

in order to put the held equations in weakly hyperbolic form. The advantage of this gauge 
over the standard Lorenz gauge is that it eliminates the derivative coupling when A = 1. 
The perturbed scalar held evolves according to the linearized equation 

□/ + Af + N :f = -4np (12) 

where 

= - (^VaV/ 3 <I> - ^F'gaf3^ , (13a) 

AT'i. = - (2AV.^$V^<I> + F") , (13b) 

play the role of external potentials and the scalar source is given by 

p = -2 j ^4(x..)rfr (14) 

:= — 2 y m(^)A(^)a(^)S 4 (x, z)dT. (15) 

The use of the “dot” and “bar” notation is borrowed from [l| and will be reviewed in a later 
section. The perturbed Einstein equation has the form 

where the perturbed Einstein tensor is given by 

2SGa,3 = -Gi'yap + 2V{a|VA7V) “ 57 '^^ - 2 .RV /3775 (16) 

+ 9(y.pF{? ' 77(5 
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The perturbed Einstein tensor is sourced by the perturbed stress-energy of the scalar field 

Svr^Ty =2V(„/V^)<h - - 9 apF'f + s<^ (17) 

- Q 7a/3 + lgaf^Fg>^^^,, 

in the bulk and the stress-energy of the point-particle 

tap = j A{^)m{^)UaUp5i{x, z)dT (18) 

on the world line. Putting everything together, the perturbed Einstein held equation takes 
the form 

□ 7 “^ -g + iV"^|./ = -167rf"^, (19) 

where 

:= 2(1 - A) - g"N ^^), (20a) 

- <5^(^V“<1>V5)$, ( 20 b) 

= -2 [(1 - X)g^^F' + 2 AV“V^$]. ( 20 c) 

Note that M“^| ^ = 0 for the choice of gauge parameter given by A = 1, which removes all 
derivative coupling from the held equations. 


III. REGULAR FIELD 

A. Field equations in compact form, Hadamard expansion, and regular/singnlar 
field decomposition 


The two coupled held equations can beiointly solved using the condensed index notation 
introduced by Poisson and Zimmerman [l|, which we briehy review here. In the condensed 
index notation, we use an upper-case latin index to denote the tensorial type of the held. 
In the present context, upper case indices will denote either a 2-tensor A = a/3 as a pair of 
indices, or a 0-tensor (scalar) A = • as the absence of tensor indices. The two helds / and 
■jap are combined into the held doublet 'i/a = {/, 7a^}- Similarly, the world line sources tap 
and p are collected into the doublet 

p^= [ g^Mix,z)q^{r) 64 ix,z)dT 
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where 


/: = 


J 4m(<l>)y4($)M“M^, A = a(3 

[ -2m(<h)v4(<l>)a(<l>), A = ■ 


( 21 ) 


and g\f{x, z) defines the parallel propagator which transports tensors at x to x' j^| 


A, A I A = aP, B'= iS' 

g B'{x,x) := { 

A = B = - 


( 22 ) 


Note that the parallel propagator is zero off the diagonal. In the condensed notation, both 
pertnrbation eqnations are combined into the single eqnation 


D^lJ^ + = -dvr/i^. 


(23) 


As fl23l) is a linear wave eqnation, its solntion can be represented as a Green fnnction 
contracted with the sonrce 

iJA = J (24) 

where Gab' obeys the wave eqnation 

aG\,{x, x') + M\V>^G%,{x, x') + N\G%{x, x') = -Ang\,5i{x, x') (25) 

and is chosen to satisfy retarded bonndary conditions, as all radiation is pnrely ontgoing. 
The principal parts of the off-diagonal Green fnnctions are solntions to homogeneons eqna¬ 
tions; e.g. DG' 1 ^^ = 0, whereas the principal parts of the diagonal Green fnnctions have 
distribntional sonrces. For all x in the convex normal neighbonrhood of a base point x', the 
retarded Green fnnction can be cast in the Hadamard form 


G^^,(x,x') = U\,{x,x')5+{a) -F I/^,(x, x')0+(-cr), (26) 

in which ct(x, x') is Synge’s biscalar, 5+(a) and 0+(—a) are the Dirac and Heaviside distribn- 
tions snpported in the fntnre of x', and U\,{x, x') and Vj^,{x, x') are smooth bitensors. The 
coefficients 17^,(x,x') and Vj^,{x,x') are recnrsively determined using the field equations 
40|. 

In the normal neighbourhood of the world line z{t) the retarded solution takes the form 
of a local leading-order r~^ piece, and a pair of tail integrals given by 

^lJ^{x) = -U'%,{x,x')q^'{u)+^ Vj^{x,z)q^{T)dT+ [ z)q^{ t) dr, (27) 

r Jr^ J-oo 
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where u is the retarded time at the point z{u) where a past directed null ray starting from 
X intersects the world line, v is the advanced time of a point connecting x to the world line 
by a future directed null ray, r is the retarded distance from x to z{u), and is the proper 


time where the world line intersects the convex normal neighbourhood of x 


40|. 


The near-zone behaviour of the solution fl27|) exhibits a Coulomb-type 1/r behaviour 
which leads to a singularity at the location of the particle r = 0. The divergence is regularized 
using the Detweiler-Whiting prescription. The method involves constructing a singular field 
V’s and subtracting it from the retarded field, i.e., V’r ■= — 'ips- The singular held is built 

from a “singular” Green function which vanishes in the causal future and past of the world 
point, is symmetric in its arguments, and solves the inhomogeneous wave equation 
Detweiler-Whiting regular held is found to be 


40|. The 


'Pr{x) = - ^^Uj^„{x,x')q^'' (v) + j Vii[x, z)q^ {t) dr (28) 

+ \ ! yii{x,z)q^{T)dT + f G\{x,z)q^{T)dT, 

^ Ju J—oo 

where radv is the advanced distance from x to the world line and v is the advanced time at 
that point. 


B. Regular field in Fermi normal coordinates 

In this section we derive expressions for the regular helds in local coordinates. We choose 
to adopt Fermi normal coordinates (FNC) adapted to the world line of the body. To con¬ 
struct a system of FNC of a point x centred on the world line, we choose a point on the world 
line X := z{t), where t is the proper time at the point. We then locate the unique spacelike 
geodesic which is orthogonal to the world line at x that connects the two points x and x in 
the normal neighbourhood. Along the connecting geodesic one can dehne a tangent bisec¬ 
tor a°‘(x,x) which is orthogonal to the four velocity m" of the world line. Here the barred 
index indicates that the quantity transforms tensorially at the barred point. The dehning 
relations for the coordinate system are therefore x^ = t, x°‘ = —e“^cT" and a°‘Ua = 0, where 
e% is a spatial triad that is Fermi-walker transported along the world line. The geodesic 
distance from x to x running along the spatial curve is given by = XaX°‘ = 2a, which is 
perturbatively small relative to the curvature scale. 

To compute the regular held in FNC, the retarded/advanced time dependencies in the 
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expression (|28]) coming from the terms x')q^' (u) and U^^„{x, x")q^'' (v) must be 

translated into dependencies on the Fermi point x. This is accomplished by starting with 
the Taylor expansion of the direct part of the Hadamard form 


U%{x,x') = g\, 


1 


<5^;, + + 0{e^) 


where 


TT^' — 
B'ti' — 


U B'fi'v' — B'u') 




(29) 

(30a) 

(30b) 


about the base point x' at the retarded time u, and then re-expanding it around the Fermi 
point X at the present time t. The retarded and advanced distances must also be expanded 
in powers of the spatial distance s with coefficients being evaluated on the world line at time 


t. The details of the procedure are spelled out explicitly in Refs. 


,Q. 


The regular field in FNC to 0{s) written in terms of condensed index notation was derived 
in |l| and found to be 


■0^(t, a;“) = —(1 — acx'^)U^{t) -|- -U^(t) dcX° -I- 0^[tail] -|- O(s^) 

o 


= -9 A 


(1 - acx‘') + g\(y^ 




BttA 


Btc 


(31a) 


(31b) 


+ 2 ^^+ V'^^hail] -t O(s^), 

where g = ^g [H, and 0^[tail] denotes the contribution from the chronological 

past 

0 "^[tail](a;) = lim f Gab{x, z{T))q^{T)dT := f Gab{x, z{T))q^{T)dT. 

J—oo J—oo 

Note that we cut the integral short before taking the limit x ^ x and not the reverse. We 
choose to work with g because it is the quantity which appears in the background equations 
of motion Amaa = gVo^* := which we interpret as the charge of the body. Under 

this definition of the charge, the “sensitivity” of the object is s := —din M/d In $ = -^Q- 
Inserting M^^ and from Sec. Ullinto Eqs. fl30l) and fl33l) we find that the U part of 
Hadamard Green function which enters the regular field is given explicitly by 

(32a) 
(32b) 


C'G = 0, 




U \.^y — -R 


6 


pv • 


(32c) 
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The gradient of the regular held, which yields the force, also includes contributions from the 
l/-terms in the Hadamard form. The Id-terms are computed with the help of the expansion 


V 


A' 

B' 


1 

4 


v" 




C /i' 


+ -N%, + 



(33) 


and we hnd 


’"“hi = - ■5“,,V-’t-Vi,* + (34a) 

= -2V“V^<h, (34b) 

^ |a/3 = “2 ~ 2^ ’ (2^^) 

1/1. = -^ (2V^4>V^4> + F") + ^R. (34d) 


The regular scalar held to next-to-leading-order in the Fermi parameter is given by the 
expression 

/r = - (1 - ttcx"") Q + (2ac + Rtc) x'' + /[tail] -h 0{s^). (35) 

To compute the regular metric perturbation we rely on the FNC expansion of the parallel 
propagator given by 


g\=l-Qax'^+ 0{s^), g^- = 0{s^), g\ = 5\ + 0{s^). (36) 

Substituting the expansion for the parallel propagator into Eq. fl31bD . we hnd that the 
regular metric perturbation takes the form 

7 r^ = - (37a) 


+ 


+ ^q^^Rtc ) + (4g*'’ac - q^^R\ct) 


X 


+ 7 “^ [tail] +0{s‘^). 


In Eqs. (15^ and (1571) . the tail terms are dehned by the relations 

7 “^[tail] := / dr z)q"^\T) + f dr G“^.(a;, 2 )g(r) (38) 

J —00 J —00 

and 

/[tail] := f drC z)q'^\T) + j drC |.(a;, 2 ;)g(r). (39) 

J —00 J —CO 
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The equation of motion is written in terms of the non-trace-reversed metric perturbation 
hai 3 = 7 a /3 “ \gai 3 l-i which we now compute. Using the expansions of the background metric 
in FNC gtt = -1 - 2aaX'^ O(s^), Qta = 0{s^), gab = Sab + O(s^), and the trace of 
7 = — (ocg^ -I- \q^^ac + ^q^^Rtc) x‘^ + O(s^), we hnd 

hft ~ ^u[fail] -|- O(s^), (40a) 

hta = Sabq^’' + ^q^^Ratctx'' + hto[tail] -h O(s^), (40b) 

hab = ~ 2 ^abq^^ + 2 ^a.b {^^o,c + —q^^^Sc + —q*^^Rt(^ x^ -\- hafe[tail] -|- O(s^). (40c) 

Using the dehnition of the background gradient Vahg^ = d^h^s — and 

the FNC expansion of the Christoffel symbols F^^^^ = aa + O(s^), F“^j = a“ -|- O(s^), with all 
other components being higher order, we hnd 

'^ahu = ^q^^cia + ^q^^Sa + ^q^^Rta + Vahtt[tail] -h 0{s), (41a) 

/ b 12 

'^thfa = Sabq “ + q^^Cla + Vthta[tSiil] + 0(s), (41b) 

and 

Vq/r = qo-a + -^q (2cia + Rta) + Va/[tail] -|- 0(s). (42) 

b 

We get an additional local contribution from the time derivatives of the tail in the normal 
neighborhood at the current time t. Recall that 

V'j^q^dr + [ G' (43) 

J —OD 

which when differentiated gives 

Vi/[tail] = + V\.q + Vi/^^^tail] + 0(s) (44) 

after using relations fIM)) and fl33l) . We use the notation Vt/R''^[tail] to denote the part of 
the tail that’s beyond leading order in the FNC expansion, which excludes the local part 
coming from V, which is smooth, and the contribution from G at the particle, which is 
singular. Likewise, we hnd 

V,fcy|tail| = A,/" (Uai,*?-'* + K.a|4) + V,fti:'ltail] + 0(s) (45) 

= - 2gV,V.* + V,hl:'|tail] + 0(a), 

where Aag^s '■= \{ga-ygp 5 + ga&gg^ - gagg^s) reverses the trace. 
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IV. EQUATION OF MOTION 


Just as we derived the field equations from the stationarity of the action under first 
variation, we find the equation of motion in the background by varying the Einstein frame 
point particle action in the background geometry 

^pp = - J dl($)m($)dr, (46) 

which gives the equation 

(mA) a“ = - (Ama) := (47) 


where 




^:= {Sj^ + u^u^) 


projects in the directions orthogonal to Ua- The evolution of the particle’s mass is dictated 
by the differential equation 


d{mA) 


= -gn^V.4>. 


(48) 


dr 

Likewise, we find the equation of motion in the perturbed geometry by varying the first-order 
perturbed point particle action 

1 


SIp = ^ J J qf dr 

with respect to the coordinate z^. We find 


(49) 


5Sl^= dr 5 A 


1 




(oq + ( qf + -q^'^h 

Working in FNC where 




D \ 


5^ 

6u^ 

D 

dr 6u^ 

leads to the equation of motion 


= 2 (f i'y’ + irffy' 


(50) 

(51a) 

(51b) 

(51c) 

(51d) 


Amtta — — (AmQ;)Va*h + qo-af + ^q Cbahtt + + q^af + -zh aq + zq ^ahu (52) 

o o o 

111 1 
- -q^^a^hab - -^q^^hta - -^q^^Vthta - -^q^^aahu 
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in FNC. Inserting the results for the regular fields and their derivatives fHOj) and making 
the substitutions, g** = AmA, M := m($)y4(<I)), and hnally Q := q{^) we hnd that the local 
self-force is given by 


{-^da + iRat " 2V„<I>Vt<I>) + 

+ QM (llaaVt$ + 2V4Va<h) + {aQ - Ma'f (53) 


Notice the similarity with the local self-force in a scalarvac spacetime. Notable though is 
the absence of the g^Va*hVt<h term. This term in the scalarvac calculation comes from 
fnVthta through the term q^^Qa, where the scalarvac mass is defined as (Js^aiarvac = rn : = 
4(m — g<h). Here, the analogous coupling in q^^Vthta has no charge squared term due to the 
difference in the dehnition of g**. Instead, we hnd a new coupling coming from the term 
ttaQ = aa{aQ — Ma')^ in the derivative of the regular scalar held. 

This expression can be expressed in order-reduced form by employing the background 
equation of motion = Q/MVa^ and its derivative da = ^'Vt'Va^ — + 

2^Vi<hVa4’. Performing the order reduction and substituting the expression for the Ricci 
tensor in terms of the background scalar held Rta = Vtd’Vad’ + 0 {s) gives 


Ft = QwJVut + 


11 


-F —Q'^ + {2Q^ - M{Ma' - aQ)) 


+ ^M{Ma' - aQ) + {aQ - Ma')'^ 
o 


~ - IqM ) 


wJu^'V^^'Vj3A> 


(54) 


for the local self-force in covariant form. 

In addition to the local force acting at the current time, we also hnd a non-local tail part 
depending on the chronological past. The contribution from within the past light cone of 
the present is given by the tail force: 




tail 


=Qaaf[tai\] + f[tail]VaQ + QVaf [tail] + ^Mahtt[tail]Va^ + ahu[tail] 

— [tail] — Ma^hta[tail] — MVthta [tail] — -Maahtt[tail], 


( 55 ) 
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which, after converting into covariant form and performing order-reduction, reads 


FT^ = 


-aQV^<h/[tail] - (Ma' - aQ) V^<h/[tail] + gv^/[tail] - [tail] 


- QV^<f>h-^4tail]u^M^ + (^V^h^5[tail] - 2V^h^5[tail]) u'^u^ 


(56) 


The perturbed regular field also contributes to the evolution equation for the particle’s 
mass 

^ = - V„<h + 2QM v?u^ - ^Q‘^R + (^Va<f>V“<h + ^F" 


+ {aQ - a'MY{u^Vo,^Y - {aQ - Ma')(M“V„<f>)/[tail] - Qu^VJ[iai\] (57) 


as well as the evolution equation for the charge 

- g'Q"(n“V„$)2 - 2MQ' (v^Vs<^ + 

+ ^QQ'R + Q"n“V„$/[tail] + g'n"V„/[tail], (58) 

where Q' = aQ — Ma' and Q” = 2a'Q + a {aQ — Ma') — Ma". Starting with initial 
conditions g(<h(0, x)) = 0, Eq. (15811 allows for a non-zero charge at some later time. In other 
words, the background scalar field can induce a charge on the world line dynamically. This 
suggests that scalarization phenomena are present in the EMRI problem as well. 


^F'g^s^ - QQ' 


V„ -F -F' 


V. SCALAR-TENSOR SELF-FORCE IN A STATIONARY BLACK HOLE BACK¬ 
GROUND 


A. Field equations 


As mentioned in the introduction, the Kerr geometry accompanied by an external, con¬ 


stant scalar field represents the mos 


to generic scalar-tensor theories 


general stationary, axisymmetric and vacuum solution 
35|. Non-stationary solutions for the scalar field configu¬ 


ration such as those arising from scalar-tensor theories with cosmological evolution or those 
with potentials leading to non-uniform asymptotic scalar field configurations may create 
hairy configurations. Though possible, their existence is highly constrained by measure¬ 
ments of period decays in binary systems displaying the absence of dipole radiation as 


predicted by hairy configurations 


20 l |. Additionally, asymptotic flatness requires that the 
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scalar potential and its first derivative with respect to the background scalar held vanish far 
from the black hole. To be consistent with a stationary black hole conhguration, we also 
demand that the background scalar held be hxed at the constant value $0 . Consequently, 
the background Ricci tensor vanishes, which implies that the stress energy of the scalar held 
is also zero. The background Einstein held equations are therefore given by the vacuum 
held equation = 0 giving rise to the Kerr family of black hole solutions for stationary, 
axisymmetric backgrounds. 

As before, we consider a non-spinning massive point-like object moving in the black- 
hole spacetime of scalar-tensor theory. The particle’s stress-energy generates perturbations 
7 q ,/3 of the background geometry. The Lorenz gauge = 0 held equations governing 

perturbations of the black hole with a constant scalar held background take the completely 
decoupled form 


□7„^ -F = -levrM J z)uaUpdT, 

□/—/i^/= —SttQ f 6i{x,z)dT, 


(59a) 

(59b) 


where M := m($o)^(‘^’o) and Q := —a(<ho)Af are the constant mass and charge of the 
point-like object, respectively, and /i^ = F"{^q) is the mass of the scalar held. The equation 
governing the metric perturbations is the well-known Lorenz gauge wave equation and the 
scalar held equation is the curved spacetime massive wave equation. 

The retarded solutions to the above decoupled system of equations are given in terms of 
the diagonal Green functions by 


70 / 3 ( 2 :) = 4M j Gaf}^ 5 {x,z{T))v?u^ dr, 
fix) = Q [ G{x,z{T))dT. 


(60a) 


Each solution is then separately regularized according to its own DW singular held to obtain 
the regular helds responsible for the self-force. The expressions for the singular/regular 
decoupled helds are widely available in the literature |40|. 


B. Equation of motion 

For a stationary black hole background spacetime with a constant scalar held = 0 
the leading order motion of the particle is geodesic in the background geometry a“ = 0. The 
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next-to-leading order motion is given by the first order self-force which is simply the sum of 
the vacuum gravitational and scalar tail self-forces 

^self = ^g“avity[tail] + ^scalar [tail] 

= (^V^h^ 5 [tail] - 2 Vyh^ 5 [tail]) u'^u^ + f [tail]. (61) 

The lack of local terms in the equation of motion is due to the geodesic leading-order motion 
of the particle and the absence of coupling terms in the linearized held equations. The mass 
evolution is dictated by the equation 

- QM"V«/[tail], (62) 

or 2 

and the charge is found to evolve according to the equation 

^ + gvv«/[tail]. (63) 

or 2 

Interestingly, we hud that the parameters characterizing particle experience a local self-held 
correction, which depends on the mass of the scalar held /i, the charge Q = and 

the prohle of the charge Q' = —M"(<I>o) at the current time, as well as a history dependent 
correction coming from the scalar held sourced by the particle in its past light cone. However, 
for realistic cosmological values, which can be as low as the Hubble scale /r 10 ^^eV, the 
local change in the mass and charge of the particle will be suppressed relative to the tail 
piece. 

VI. CONCLUSION AND FUTURE DIRECTIONS 

In this paper we derived the regular held and the equation of motion for a massive point- 
particle in the Einstein frame of scalar-tensor theory to hrst-order in the mass-ratio. We hnd 
that the self-forced evolution in scalar tensor theory is diherent from the evolution of a mas¬ 
sive, scalar charge in an scalarvac spacetime of general relativity. For non-stationary scalar 
held backgrounds, the leading-order motion is accelerated in the background spacetime and 
the hrst-order self-force correction takes the form of a set of local and non-local terms. The 
non-local, “tail”, component of the force is characterized by non-diagonal Green functions 
which result from the coupling in the held equations. The local ST self-force resembles 
the local self-force on a charged, massive particle in a scalarvac spacetime with additional 
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couplings due to the time variability of the particle’s charge in addition to its time variable 
mass. The reason is that in ST theory, the strong equivalence principle is violated and the 
motion is sensitive to the internal constitution of the body. Here we used the point-particle 
approximation, where the details of the body’s internal structure are coarse grained out. 
However, the mass and charge of the object carry a scalar held dependence, and this creates 
fundamentally different world line dynamics. Furthermore, the evolution of the mass and 
charge of the particle is inhuenced both locally and non-locally by the regular self-held. This 
is true even for black holes with stationary scalar backgrounds, where the mass and charge 
of the particle evolve locally if the scalar held has mass and non-locally through the regular 
part of the scalar perturbation. This has consequences for the boating orbit scenario. It 
is evident from our equations that self-held corrections to the particle’s mass may cause a 
violation of the super-radiance condition or lead to an imbalance in the radiation huxes, 
precluding any prolonged dynamical boating of the orbit. Cardoso et ah estimated the 
ehect of ingoing negative scalar radiation hux at on the mass of the large black hole, 

but they neglected to mention the possibility for the scalar particle to possess a dynamical 
mass. To estimate whether the particle’s mass variation is important we can compare the 
timescale on which the particle’s mass varies to the timescale on which the black hole’s mass 
secularly decreases. When the orbit is in a boating conhguration, the negative scalar hux at 
the event horizon is balanced by the hux of outgoing gravity waves at ■ This cannot last 
forever due to the decrease in the black hole’s mass as the scalar held extracts energy from 
it. The timescale for the extraction is roughly Tbh ~ where Mbh is the mass of 

the large Kerr black hole and M is the mass of the particle. The timescale governing mass 
changes to the small body due to the self-force is roughly Tgb = M/(Q^/i^) ~ l/(M/i)^ when 
ignoring the tail contribution. For self-force ehects to be important they must occur before 
or at least around the time that signihcant mass loss of the black hole from the negative 
hux occurs, or Tgb ^ Tbh or M/Mbh ^ Thus, self-force ehects become important for 

a large region of the parameter space. 

The self-force for ST black holes with trivial scalar conhgurations is also diherent than in 
GR. In vacuum GR, the force is completely non-local and involves the tail force due to the 
metric perturbation only. In scalar-tensor theory, even with a constant background scalar 
held solution, the scalar perturbation will also inhuence the motion: a point particle orbiting 
an ST black hole experiences both the metric and scalar tail forces. These deviations from 
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the GR self-force will likely have observational consequences. The self-force is known to 
cause a shift in the innermost stable circular orbit, alter the rate of periastron advance 
for eccentric orbits and influence the renormalized redshift of photons sent from the 


particle to distant observers 


m 


Next steps include computing observational effects in ST theory and using the results to 
constrain its parameters. For instance, one can imagine setting up a numerical calculation 
of the self-force on a point-particle around a star or black hole in ST theory using the mode- 
sum method. This will require formulating the coupled held equations in a fashion where the 
modes of the retarded held can be numerically determined and the mode representation for 
the singular held can be subtracted. Linz and collaborators have shown that the derivatives 
of the singular helds for a scalarvac spacetime completely decouple and that the regulariza¬ 


tion parameters for the force are the sum of the decoupled parameters 


32| . It is reasonable 


to expect that the singular held gradients will also be decoupled in the present case, as the 
equations have the same form, and this will certainly alleviate some of the computational 
challenge as only the retarded held will need to be solved from coupled equations. One 
could also imagine solving the decoupled held equations fl5^ to compute the self-force and 
mass/charge evolution for a particle outside a Kerr black hole in ST theory. One could then 
choose initial data and scalar held parameters corresponding to a boating orbit of Ref. 271 
and compute the self-force and mass evolution to assess the stability of the orbit. 
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